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Abstract 

Using the isospin dependent quantum molecular dynamics model, we study the effect of charge 
asymmetry and isospin dependent cross-section on different aspects of elliptical flow. Simulations 
have been carried out for the reactions of ^'^'^Xm -|-^^^ Xm, where m = (47, 50, 53, 57 and 59) and 
'^^Xn+'^^ Xn, where n= (14, 16, 18, 21 and 23). Our study shows that elliptical flow depend strongly 
on the isospin of cross-section. The transition energy remains almost constant with increase in N/Z 
of the system. A good agreement is obtained with experimental measurements. 
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I. INTRODUCTION 



Heavy-ion collisions provide a possibility to study the properties of nuclear matter in 
conditions that are vastly different than reported in normal nuclei l| . Considerable progress 



has been made recently in determining the equation of state of nuclear matter from heavy ion 
reaction data [2,]. A prominent role among available observables is played by the collective 



flow. Much theoretical and experimental efforts have been made to the study of collective 
flow in HICs 3|. The elliptic flow has proven to be one of the most fruitful probes for 



extracting the EoS and the dynamics of heavy-ion co^ 



lisions. The parameter of elliptic flow 



is quantified by the second-order Fourier coefficient ^ from the azimuthal distribution of 



detected particles at midrapidity as: 



^ = + 2viCos(f) + 2v2Cos2(j)). 



(1) 



Where is the azimuthal angle of emitted particle momentum relative to the x axis. 
Positive values for ( Cos2(j) ) refiect a preferential in-plane emission, whereas negative values 



of ( Cos2(l) ) refiect a preferential out-of-plane emission. 



Yu-Ming Zheng et al., 



5| studied the proton elliptic fiow in the collisions of ^^Ca Ca 



at energies from 30 to 100 MeV/nucleon. They showed that, with increasing incident energy, 
the elliptic fiow shows a transition from positive to negative values. Its magnitude was 



found to depend on both nuclear equation of state (EOS) as well as on the nucleon-nucleon 
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scattering cross-section. On the other hand, D. Persram et al., 6| have shown, for the first 



time, the effect of parameter-free self-consistent calculation of nucleon-nucleon in-medium 



scattering cross-section implemented BUU model on different flows. Their results favored 
in- medium cross-section compared to free one. Y. Zhang et al., |7| have investigated 



the elliptic flow in heavy-ion collisions at energies from several tens to several hundreds 
of MeV/nucleon. They showed that, soft nuclear equation of state and incident energy 
dependent in-medium nucleon-nucleon cross-sections are required to describe the excitation 
function of the elliptic flow at intermediate energies. 



II. ISOSPIN-DEPENDENT QUANTUM MOLECULAR DYNAMICS (IQMD) 
MODEL 



Our study is performed within the framework of IQMD lOj model where hadrons prop 



agate with Hamilton equations of motion: 



drj 
dt 



d{ H ) 
dpi 



dpi 
dt 



d{ H ) 
dri 



(2) 



with 



H) = {T)+{V) 



i * j j>i 

X fj {r, p, t)drdr' dpdf) . 



(3) 



3 



The baryon-baryon potential V^^ , in the above relation, reads as: 



V'^^(r' — r\ = V^^ 4- V'^ -+■ V'^ + V"-^ 
V' '/ Skyrme ' Yukawa ' Coul ' sym 



+ 13 — 71—; / X — + 



(IP — r|//x) 1^ — r| 

^t,-TlTib{r{-r-). (4) 

Here Zi and Zj denote the charges of i*^ and j*'' baryon, and T3, T| are their respective 
T3 components (i.e. 1/2 for protons and -1/2 for neutrons). Meson potential consists of 
Coulomb interaction only. The binary nucleon-nucleon collisions are included by employing 
the collision term of well known VUU-BUU equation. During the propagation, two nucleons 
are supposed to suffer a binary collision if the distance between their centroids 



\ri-rj\< (^tot = (^(Vs, type) , (5) 

"type" denotes the ingoing collision partners (N-N, N-A, N-tt,..). In addition, Pauli blocking 
(of the final state) of baryons is taken into account by checking the phase space densities in 
the final states. 



III. RESULTS AND DISCUSSION 

For the present analysis, simulations are carried out for two sets of the reaction using 
soft equation of state. For the first case, the mass of the colliding nuclei is fixed to be 124 
units and for the second set, the mass of the colliding nuclei is fixed to be 40 units. In other 
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words, we study, the reactions of ^^^X^ +^^^ X^, where ^^"^X^ = {^'^'^Ag^r, ^"^^Sn^o, ^^^/sa, 
^^^La57 and ^^^Prsg), respectively. The second set corresponds to: H-^" y„, where ^''Y'^ 
= (^°^3) '^''5'c2i, "^"Aris, "^^Sie and ^^Sz^), respectively. The phase space generated by the 



10|] method. The 



IQMD model has been analyzed using the minimum spanning tree (MST) 
elliptical flow is defined as the average difference between the square of x and y components 
of the particles transverse momentum. Mathematically, it can be written as: 

2 2 

{v,} =< Cos2ct^ >= (%^) (6) 

where and Py are the x and y components of the momentum. The positive value of 
elliptical flow describes the eccentricity of an ellipse-like distribution and indicates in-plane 
enhancement of the particle emission. On the other hand, a negative value of f 2 shows the 
squeeze-out effects perpendicular to the reaction plane. Obviously, zero value corresponds to 
an isotropic distribution. Generally, for a meaningful understanding ( f 2 ) is extracted from 
the mid-rapidity region only. Naturally, mid-rapidity region (-0.1 < X"'"" < 0.1) corresponds 

^beam- 

to the collision (participant) zone and hence signifies compressed matter. On the other hand, 
" 7^ corresponds to the spectator zone, ( J"-"^ ( -0.1) corresponds to target like (TL) 

matter and ( J""-" ) 0.1) corresponds to projectile like (PL) matter. 

To study the effect of charge asymmetry on the elliptical flow as a function of { Pt ), 



(where ( ) is transverse momentum of particle and is given by pt = \J {pi +Py))- We 
display, final state elliptical flow in Fig.l and 2 for free particles (A = 1) (upper panel), 
IMF's (2 < A < 4) (middle) and IMF's (5 < A < Atot/Q) (lower panel) at an incident 
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FIG. 1: Transverse momentum dependence of the elliptical flow, summed over the entire rapidity 
distribution, at 6 = 0.5 for three different reactions at 50 (left) and 100 (right) MeV/nucleon. 

energy E = 50 MeV/nucleon (left) and E = 100 MeV/nucleon (right) for the reactions of 
+124 where m = 47 and 59 (in Fig.l) and 4°F„ +4° where n = 14 and 23 (in 
Fig. 2). Here elliptical flow is summed over all rapidity bins. Figs.l and 2 reveal: 
(a) Gaussian shape is obtained for ( ) in all cases. Note that the elliptical flow is integrated 
over the entire rapidity range. This Gaussian shaped behavior is quite similar to the one 



reported by Colona and Di Toro et al., |ll |. 



(b) Peak of the Gaussian shifts towards lower value of ( Pj ) for the heavier fragments. This 
is because the free nucleons and LMF's feel the mean field directly, while heavy fragments 
have weaker sensitivity 12 1. 



c) The neutron rich systems Ag4i+ Ag^-j and Siu+ Siu exhibit weaker squeeze-out 
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FIG. 2: Transverse momentum dependence of the elliptical flow, summed over the entire rapidity 
distribution, at 6 = 0.5 for three different reactions at 50 (left) and 100 (right) MeV/nucleon. 

flow compared to the neutron deflcient reactions ^"^^Pr^g +^^^ Pr^g and ^^^¥23 +'^^ V23. Our 



findings are in agreement with Zhang et al., 13|] where a neutron rich system was found to 
exhibits weaker squeeze-out flow. 

To study the effect of isospin dependence of cross-section and charge asymmetry on the 
elliptical flow, we display in Figs. 3-6, the transverse momentum dependence of elliptical flow 
for the reactions of ^^^Ag^j +^^4 ^^^^^ ^^^Pr^g +^^^ Pr^g, ^^Siu ^^i4 and ^°V^23 V23. 
We divided total elliptical flow into contributions from target-like (TL) (left panels), 
mid-rapidity (middle panels) and projectile-like (PL) (right panels) particles at E = 100 
MeV/nucleon. The upper, middle and bottom panels represent the free nucleons, LMF's 



and IMF's. These flgures reveal following points: 
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FIG. 3: Transverse momentum dependence of the elliptical flow at E = 100 MeV/nucleon for two 
different cross-sections for the reaction ^'^'^Ag^j +1^4 j[g^.^_ 
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FIG. 4: Same as Fig. 3 but for the reaction ^'^^Pr^g +^'^^ Pr^g. 
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FIG. 6: Same as Fig. 3 but for the reaction V23 + V23. 
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FIG. 7: Variation of elliptical flow with energy for the reactions ^^^X^n +^^^ where ^^^X^ = 
(124^547, '''Sn,o, '''h3, '''La,, and ''^Pr,,). 

(a) ( f 2 ) is sensitive to different nucleon-nucleon cross-sections. Weaker squeeze-out ffow 
is observed in the case of isospin independent cross-section. This happens because in the 
case of isospin dependent cross-section, neutron-proton cross-section is three times larger 
compared to neutron-neutron and proton-proton cross-section that will enhance binary 



collisions. These findings are in agreement with ref. 14 1. 

(b) Comparison of Figs. 3 & 4 and Figs. 5 & 6 indicates significant dependence of ( f 2 ) on 
;he charge asymmetry of colliding pairs. Our findings are also supported by Zhang et al., 



13| 



(c) Closer looks to these figures help us to understand the origin of these isospin effects. 
From the figures, it is clear that percentage change is maximum in the case of mid-rapidity 
region for both isospin dependent and isospin independent cross-section. Indicating that 
isospin effects originates from mid-rapidity region. 



In figs. 7 and 8, we display the variation of excitation function of elliptical flow ( V2 ) 
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FIG. 8: Variation of elliptical flow with energy for the reactions "^^Yn +'^^ Yn, where '^^Yn = (^''V23, 
40^021, ^°Ans, '°5i6 and ^^Siu). 

for free nucleons, LMF's and IMF's for mid-rapidity region using same set of the reactions 
considered earlier. We note: 



(a) The elliptical flow turns negative with 
faster after ( f 2 ) reaches a minimum value 



3eam energy. This is because spectators move 



15|. The energy at which this behavior changes 



is found to decrease with the size of the fragment. This means that the flow of heavier 



fragments is larger than that of LM 



findings are in agreement with ref. 16 1 



"s and free nucleons at all beam energies. These 



(b) There occurs a transition from in-plane emission to out-of-plane. This is due to the 
fact that the contribution of participant zone dominates the reaction in midrapidity region 
leading to the transition from in-plane to out-of-plane. In other words, participant zone is 
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FIG. 9: Transition energies {ETrans in MeV/nucleon) for elliptical flow as a function of N/Z. 
primarily responsible for the transition from in-plane to out-of-plane. The energy at which 



this transition observed is dubbed as transition energy ETrans- 



Now to study the effect of charge asymmetry on the transition energy of free nucleons, 
LMF's and IMF's, we show in Fig.9, the transition energy as a function of N/Z for two 
sets of reaction. From both the figures, it is clear that transition energy remains almost 
constant with increase in the N/Z of the system. 



To further strengthen our interpretation of the results of elliptical flow ( f2 ), we display 



in Fig. 10 a comparison of theoretical results of elliptica 



tracted by the INDRA@(GSI+GANIL) collaboration 



flow with experimental data ex- 



19]. Here simulations are performed 
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FIG. 10: The energy dependence of elliptical flow for the reaction ^^^Snso +^'^^ Sn^Q. 



124 



for the reaction 5*^50 + Sn^Q with cTj^o reduced by 20%. It is worth mentioning that 
the results with above choice of cross-section are in good agreement with the experimental 
data of ref. 17|. The choice of reduced cross-section has also been motivated by ref. ij] 



as well as many previous studies 



19|. 



IV. SUMMARY 

Using the isospin dependent quantum molecular dynamics model, we have studied the 
effect of charge asymmetry and isospin dependent cross-section on different aspects of 
elliptical flow. Here simulations have been carried out for ^^"^X^ +^^'^ X^, where m = 
(47, 50, 53, 57 and 59) and ^0X„ X„, where n= (14, 16, 18, 21 and 23) . Our study 
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shows that eUiptical flow is observed to be strongly dependent on the isospin-dependent 
cross-section. The transition energy remains almost constant with increase in N/Z of the 
system. The comparison with experimental data of INDRA@(GSI+GANIL) collaboration 
supports our findings. Moreover, our results of will be of great use for the 

experimentahsts working at SCC500 at VECC Kolkata (INDIA). 
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